Many-body Landau-Zener Transition in Cold Atom Double Well Optical Lattices 
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Ultra-cold atoms in optical lattices provide an ideal platform for exploring many-body physics of 
a large system arising from the coupling among a series of small identical systems whose few-body 
dynamics are exactly solvable. Using Landau-Zener (LZ) transition of bosonic atoms in double well 
optical lattices as an experimentally realizable model, we investigate such few to many body route by 
exploring the relation and difference between the small few-body (in one double well) and the large 
many-body (in double well lattice) non-equilibrium dynamics of cold atoms in optical lattices. We 
find the many-body coupling between double wells greatly enhances the LZ transition probability, 
while keeping the main features of the few-body dynamics. Various experimental signatures of 
the many-body LZ transition, including atom density, momentum distribution, and density-density 
correlation, are obtained. 

PACS numbers: 03.75.Lm, 05.70.Ln 



Understanding many-body physics in strongly- 
correlated lattice models is essential for the explanation 
of many important condensed matter phenomena, such 
as the high temperature cuprate superconductivity [![. 
In this context, ultra-cold atoms in optical lattices 
provide an ideal platform for emulating numerous phe- 
nomena in solids because of their ability of accurately 
implementing various lattice models without impurities, 
lattice phonons, and other complications 2-6]. In solids, 
a large many-body system may be composed of a series 
of small identical few-body systems, therefore it would 
be natural and interesting to investigate how many-body 
properties (e.g. correlations) of the large system emerge 
or differ from the few-body properties of the small 
systems (see Fig. [TJ, for an illustration) Q. Although 
such few to many body route may provide a unique angle 
for understanding the underlying many-body physics, 
its experimental realization is very challenging in solids. 
In contrast, such route may be easily explored using the 
recent experimentally realized cold atom double well 
optical lattices [8l4l0|. where the few-body dynamics 
in each double well can be solved exactly, while the 
many-body physics emerges from the inter-well coupling. 

The double well lattices not only allow studying in- 
teresting many-body ground states 111- 1A \ , but also the 
equilibrium and non-equilibrium dynamics of cold atoms 
after an adiabatic or sudden change of the atom or lattice 
parameters [lH, [l|| . Recently, the non-equilibrium dy- 
namics in optical lattices after a sudden quench have been 
investigated intensively IB"22|. While the dynamics for 
the adiabatic process are expected to follow the change of 
the system Hamiltonian, the goal for studying the quench 
dynamics is to understand the non-equilibrium physics 
and the relaxation to the equilibrium steady states in 
the presence of many-body interactions. Although the 
two limiting cases (adiabatic or sudden) have been widely 
studied, the intermediate region, that is, the parameter 
variation with a finite rate, has been largely unexplored. 

In this Letter, we integrate these two important as- 
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FIG. 1: (a) Illustration of the route from the few-body physics 
in identical small systems to the many-body physics in a large 
system through coupling, (b) An example for (a): the LZ 
transition of cold atoms in ID double well optical lattices. The 
coupling between different double wells (the small systems) is 
provided by the inter- well coupling J'. 



pects (i.e., route from few to many body physics and non- 
equilibrium dynamics) for cold atom optical lattices into 
one simple, but experimentally feasible model: the many- 
body Landau-Zener (LZ) transition in a one-dimensional 
(ID) double well optical lattice. In the LZ transition, 
the parameters of the Hamiltonian vary with a finite 
rate (neither adiabatic nor sudden), and the dynamics 
are naturally non-equilibrium. Furthermore, the dynam- 
ics of a single atom in an isolated double well, a classical 
example of the LZ transition (23[, are exactly solvable. 
Recently, the LZ transition has been generalized (both 
theoretically and experimentally) to a BEC in a double 
well with the mean- field nonlinear interaction (25l - |28| . 

Here we investigate the many-body LZ transition in a 
double well optical lattice to explore the route from few 
to many body non-equilibrium dynamics. We mainly fo- 
cus on the following issues: (i) the LZ transition for a 
few interacting atoms (in contrast to hundreds of atoms 
in the mean- field region |25H27| ) in an isolated double 
well (i.e., without inter- well coupling); (ii) the collective 
LZ transition of many atoms in the double well lattice 
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with the inter-well coupling. Our goal is to explore the 
difference as well as the relation between the few-body 
physics in (i) and the many-body physics in (ii) . We find 
that the onsite interaction U (repulsive or attractive) be- 
tween atoms can strongly modify the LZ transition in 
(i). For the repulsive interaction, there is an oscillation 
of the LZ probability with respect to U. The inter-well 
tunneling that couples different double wells can signifi- 
cantly increase the LZ transition probability, while still 
keep the main features of the LZ transition process in a 
single double well. We show the signature of the many- 
body LZ transition in various experimentally measurable 
quantities such as the atom density and momentum dis- 
tributions, and the density-density correlation. 

The Hamiltonian: The ID double well lattice, 
schematically shown in Fig. [TJ), has been realized in 
experiments by superimposing two optical lattices with 
different wavelengths [8l4ld|. Within the tight-binding 
approximation, the dynamics of atoms in the ID lattice 
are described by the Bose-Hubbard model [5[, 
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where a, (aj) are the annihilation (creation) operators 
for bosons at the lattice site i, nt = ajdj, J (</') repre- 
sents the intra-well (inter-well) tunneling, 2A = 2-ft is 
the chemical potential difference between two neighbor- 
ing lattice sites in a single double well, which varies with 
the time t with the sweeping rate 7. U is the on-site 
interaction strength between atoms. 

The time-dependent dynamics of atoms governed by 
Hamiltonian ([1} is solved numerically. Because the di- 
mensionality of the Hamiltonian increases exponentially 
with respect to the lattice size, we simulate the dynamics 
of this system using the recently developed time-evolving 



block decimation (TEBD) algorithm [29H32j | with an open 
boundary condition, which is a powerful tool for study- 
ing lattice dynamics in one dimension. For a small lattice 
size, we also study the dynamics using the Runge-Kutta 
method and find excellent agreement with the TEBD 
method. Henceforth, the energy unit is chosen as the 
intra-well tunneling J, and the corresponding time unit 

is hj- 1 . 

Few-body dynamics in one double well: For cold atoms 
in isolated double wells without inter- well tunneling (i.e., 
J' = 0), the probability for the atoms remaining on the 
same state after they pass through the LZ transition 
regime is known to be exp (— ttJ 2 /^ 23[ for U = 0. In 
the presence of the on-site interaction U for several atoms 
in one double well, the LZ transition can be dramatically 
different because the interaction shifts the energy lev- 
els for atoms. For simplicity, we illustrate the essential 
physics using two atoms (N = 2) in one double well. Un- 
der the number state basis {|2l0_r) , |1l1_r) , \0l2r)}, the 
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FIG. 2: (Color online) Plot of the instant energy levels in a 
single double well with two atoms. The level crossing depends 
strongly on the sign of U. The arrow indicates the time evo- 
lution of the quantum state for a relatively small sweeping 
rate 7. 



Hamiltonian can be written as 



Ho 




(2) 



where L (R) represents the left (right) well. Note that 
the quantum dynamics of the above simple Hamiltonian 
cannot be solved analytically using the S'-matrix method 

The instant eigenenergy levels at the time t are plot- 
ted in Fig. O With the large \U\ > J, the two large 
anticrossings at t 7^ with the gap ~ J yield direct LZ 
transitions between two quantum states that differ by 
exactly one atom, while the tiny anticrossing at t = 
with the gap ~ J 2 j\J corresponds to the indirect second 
order transition process, which has a negligible contri- 
bution to the total LZ transition probability when the 
sweeping rate is faster than the time scale determined by 
this gap. Assuming initially the atoms are prepared on 
the ground state \2l0r) for t —> —00, we see for U > 0, 
all atoms can be transferred from the left to the right well 
at t —> +00 when 7 is small (the arrow in Fig. 2a). While 
for U <C — J at most one atom can be transferred for a 
relatively slow sweeping rate 7 (but still fast enough such 
that the tiny gap at t = can be neglected). For a very 
small sweeping rate 7 (thus the gap at t — cannot be 
neglected), the dynamics are still adiabatic and all atoms 
can be transferred to the right well, as expected. For the 
atom number larger than two (N > 2), the physical pic- 
ture shown in Fig. [2]is still similar. Generally, for U 3> J 
we observe N different LZ transitions, while for U <C — J 
only one direct LZ transition can be found (hence at most 
one particle can be transferred from left to right if the 
sweeping rate is faster than the indirect transition gap). 
Finally the time interval between two adjacent direct LZ 
transitions is exactly the same St = 2U/'y. 

Many-body dynamics in double well lattices: Now we 
switch on the coupling J' between double wells to in- 
vestigate how J' and U influence the LZ transitions. To 
illustrate the essential physics in the large lattice, we first 
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FIG. 3: (Color online) Plot of the atom number 113 in the 
third site versus (a) J' and (b) U for two coupled double 
wells. 7 = 2ty. 



consider an isolated system with only two coupled double 
wells, each of which contains N — 2 atoms. Similar as 
the Hamiltonian @ for a single double well, we rewrite 
the Hamiltonian ([TJ to a 35 x 35 matrix form under the 
number state basis {\i1i2i3H)}, where ii — {0,1,2,3,4} 
represents the number of atoms at the site I with the con- 
straint i\ + 12 + 73 + 74 = 4. We then numerically solve the 
dynamics of atoms using the fourth-order Runge-Kutta 
method to obtain the exact dynamics. The quantum 
state of the system is initially prepared at \1l^>r)®\^l^r) 
(the ground state at t — > —00), and the number of atoms 
in the third site 773 for a large positive t (effectively 
t — > +00), i.e., the number of atoms remaining in the 
left well, is calculated. Note that a small 713 means a 
large LZ transition probability. 

In Fig. [3^, we plot 773 as a function of J' with the 
finite U. 713 decreases monotonically with the increas- 
ing J' for all interaction regions, which indicates the LZ 
transition probability is enhanced by J'. This can be un- 
derstood from the fact that when the inter- well tunneling 
is switched on, the atoms in the third site can tunnel to 
not only the fourth site, but also the second site. 

In Fig. [3j3, we plot 773 with respect to U for finite 
J'. We see a strong dependence of 713 on the sign of U. 
For a large attractive U, 773 is large and the LZ transition 
probability is small because atoms are bound together for 
tunneling by the large attractive interaction, as shown in 
Fig. 2b. The inter-well tunneling J' barely modifies the 
transition probability. While in the repulsive U region, 
773 decreases quickly with the increasing U, but saturates 
at the large U limit. In the large U region, 773 oscillates 
periodically with U even for J' — 0. Note that there are 
two direct LZ transitions for two atoms in a double well. 
After passing the minimum gap point of the first LZ tran- 
sition, 713 still oscillates periodically with time. Since the 
time interval between two adjacent direct LZ transitions 
is St = 2J7/7, it is expected that 713 before the second LZ 
transition may oscillate with U, leading to the final oscil- 
lation dependence of 713 on U at t — > +00. The nonzero 
J' reduces 773 significantly by enhancing the LZ transi- 
tion probability. However, the positions of the peaks and 
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FIG. 4: (Color online) (a) Plot of the atom density distribu- 
tion m for differentlattice sizes M. J' — 0.5, U = 5, 7 = 2tt. 
(b) Plot of the number of atoms at the left site ul of the 
central double well of the lattice with respect to U. M — 25, 
and 7 = 2-7T. 



valleys of the oscillation are the same for different J'. 
Note that general analytic perturbation methods (with 
J' as a small parameter) do not work for the coupled LZ 
transition because of the non-adiabaticity and the time- 
dependence of the system, where even the dynamics in 
one single double well cannot be obtain analytically. 

With the knowledge of the LZ transition in single and 
coupled double wells, we now study the many-body LZ 
transitions in a lattice with many (~ 30) double wells 
using the TEBD algorithm. Note that because of the 
exponentially large dimensionality of the lattice system, 
the fourth-order Runge-Kutta method (as well as any 
exact diagonalization method) does not work. The initial 
wavefunction is assumed to be \0 l r )® Nl ® \2 L R )^ M <g) 
IOlOa)®^, where M is the number of occupied double 
wells. Nt {Nr) is the number of unoccupied double wells 
at the left (right) edge of the lattice to eliminate the 
effects of the open boundary condition. With the finite 
J', the atoms may diffuse to the left and right edges in the 
time scale ~ 1/J'. For a relatively large 7, such diffusion 
does not affect the LZ transition at the center of the 
lattice, as we show below. In the numerical simulation, 
we choose Nl = Nr = 3, which are sufficient for a wide 
range of 7. 

To check the effects of the atom diffusion on the LZ 
transition, we plot the number of atoms (at t — > +00) in 
each lattice site in Fig. 0k for 7 = 27r and different M. 
Although the atom diffusion at the edge of the lattice is 
significant, the atom number fluctuation in the central 
part of the lattice is generally very small. For instance, 
the difference between the final atom number at the cen- 
ter of the lattice is smaller than 0.002 for two different 
lattice sizes M = 23 and 31. Our numerical results show 
that when M > 20 the LZ transition in the central part 
of the lattice can be a good approximation for that in an 
infinite large lattice M — > 00. 

We calculate the number of atoms rij, (at t — > +00) 
in the left site of the double well at the center of the 
lattice with M = 25, and plot it as a function of U 
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FIG. 5: (Color online). Plot of the atom density rn{t) (al, a2), 
momentum distribution n(k,t) (bl, b2), and density- density 
correlation s(k,t) (cl, c2). J' — 0.5. The left column: U = 
— 5; the right column: (7 = 5. 7 = 2n. 



in Fig. |4}d. To compare the results with the few-body 
physics in the small system, we also present the results 
for single and two coupled double wells in the same figure. 
When U < 0, nj, is generally larger than one because of 
the reduced LZ transition probability by the attractive 
interaction, is almost the same as that in coupled 
double wells, which indicates J' does not modify the LZ 
transition significantly for the attractive interaction. For 
U > 0, riL is generally smaller than 1. Although in a 
large quantum system the instant eigenenergy levels be- 
come very complex, we still observe the oscillation of ul 
with respect to U . However, the perfect oscillation in 
the coupled double wells is smeared by the interference 
effects in lattices, leaving only one big dip at U ~ 4.5. 

Experimental observation: The many-body LZ tran- 
sition can be observed using several different methods, 
(i) The number of atoms in the left or right wells of 
the lattice (e.g., in Fig. 0} can be measured us- 
ing the band mapping method (loj . In Figs. [5^, we 
plot the atom density distribution during the LZ tran- 
sition in the central part of the lattice. We see there 
is only one LZ transition for U -C — J, while several 
LZ transitions for U S> J. (ii) When one atom is 
transferred from one site to the neighboring site, there 
should be a change of the atom momentum distribu- 
tion n(k,t) = J2i j e lki - l ~ti {ip\a\a 3 ;\ip), which is plotted 



in Fig. [5}d and can be measured using the time-of-flight 
image. Here ip (t) is the wavefunction obtained from the 
TEBD numerical simulation. Note that the actual mo- 
mentum distribution may oscillate rapidly after the LZ 
transition, which may not be observable in experiments. 
Here, we take the average momentum distribution in a 
small time interval 0.9S.J" 1 in Fig. [SJd. Clearly, the atom 
momentum distribution changes dramatically during the 
LZ transition. There is a one-to-one correspondence be- 
tween the density and momentum changes in the LZ 
transition, (hi) The density-density correlation s(k,t) = 
e lk{l ~ j) [(ip\riinj\ij}) - (MnM)(MnM)], which can 
measured using the Bragg scattering [33], shows similar 
features because niUj should also have a sudden change 
in the LZ transition process. Note that n(k,t) shows a 
strong asymmetry about k — because the momentum is 
transferred along a specific direction during the LZ tran- 
sition. In contrast, the density-density correlation, due 
to its definition, is exactly symmetric about k = 0. 

Experimentally, the ID double well lattice has been re- 
alized for 87 Rb atoms by superimposing two optical lat- 
tices with the wavelength Ai = A2/2 and the correspond- 
ing lattice potential depths V\ and V2 [l(|. The motion 
of atoms along the other two dimensions are frozen by 
two additional regular lattices with the potential depth 
V± = ZQE R , where E R = h 2 /2m\ 2 = 2tt x 3.7 KHz is 
the atom recoil energy. By tuning V\ and V2, we can set 
the intra- well tunneling energy J ~ h x 2ir x 200 Hz. The 
typical onsite interaction strength in the optical lattice 
is U ~ h x 2-7T x 1 KHz, but can be tuned using Fes- 
hbach resonance [3J|. The chemical potential A can be 
tuned by shifting one laser beam with respect to the other 
[Ioj |. The typical LZ transition time is ~50?i/J ~ 40 ms, 
which can be easily achieved in experiments. Taking into 
account of the symmetry of the energy levels for the re- 
pulsive and attractive interactions (Figs. 2a and 2b), we 
can use the repulsive interaction (U > 0) to emulate the 
LZ transition with the attractive interaction (U < 0) by 
initially preparing the atoms in the excited state (the up- 
per branch in Fig. 2a), similar as the method used in a 
recent experiment [Toj . 

Conclusion: In summary, we study the many-body ef- 
fects in the LZ transition using ultra-cold bosonic atoms 
in ID double well optical lattices. The dynamics of such 
non-adiabatic systems are obtained through large scale 
numerical simulations using the TEBD algorithm. Our 
results arc important for understanding not only the 
route from few to many body physics when individual 
small systems are coupled to form a large system through 
the many-body coupling, but also the non-equilibrium 
dynamics in optical lattices when the system parameters 
vary with a finite rate (neither adiabatic nor sudden). 
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